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Notations 


The following notations has been used in this. Thesis 

a,; Acceleration in x-direction 

ay Acceleration in y-direction 

Sz Acceleration in z-direction 
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[J] Jacobian matrix 

[K] Stiffness matrix 

[K®] Element stiffness matrix 
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L Length of pipe 

[L] Matrix relating strains and displacements 

Mg Confined modulus 
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Shape function matrix 
Number of elements 
Intensity of surface loading 
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Potential energy 
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Time at nth step 
Kinetic energy 
Displacement in x-direction 
Relative displacement 
Absolute displacement 
Approximating function over element 
Ground displacement 
Displacement component 
Ground acceleration 

Mode transformation 
Three dimensional relative displacement 
Three dimensional relative velocity 
Three dimensional relative acceleration 
Nodal displacement vector 
Displacement in y-direction 
Displacement in z-direction 
Nodal coordinates 
Vector of nodal coordinates 
Maximum response of ith node 
Vector of normal coordinates 
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Strain over element 
Surface of domain 
Surface of the element 
Volume of the domain 
Volume of the element 
Poisson’s ratio 

Poisson’s ratio of the pipe material 
Poisson’s ratio of the soil 
Matrix of mode shapes 
Frequency 
Mass density 
Cross-modal coefficient 
Damping ratio for mode i 



Abstract 


Pipelines as a part of lifeline systems are very much important in normal life and 
particularly so during and after earthquake. Basic services such as water supply, sewage 
disposal, gas and petrol supply, electrical and telephone communications etc. may suffer 
severe damages due to disruption caused to the ground installations including buried pipes. 
'Fhe seismic behaviour of buried pipes is quite different from that of super structures 
mainly because of inertia forces and ground motion. 

The present study is focussed on the response analysis of buried pipes under seismic 
loading using a 3-D finite element model. Finite element method package, SAP-80 has 
been used to analyze the pipe-soil system imder different load conditions and different soil 
conditions. Parametric study has been conducted on 1) response of pipe-soil system under 
static loads and 2) response of pipe-soil system under seismic loads. The parameters 
considered are - load width-to-radius ratio, radius-to-thickness ratio, embedment depth-to- 
diameter ratio, mode shape number, modular ratio, 3-D seismic loading and liquified soil 
zone. The results obtained, have been presented for a wide range of parameters and are 
compared with those already reported in literature wherever possible. 

It has been observed that the present 3-D thodei can be used effectively to obtain the 
response of pipe-soil system undeir dilfetent k))Sds. 



Chapter 1 


Introduction 

1.1 General 

Safe and reliable functioning of buried pipelines , as a part of lifeline systems are 
very much important in normal circumstance and particularly so during and after 
earthquake. For instance, basic services such as water supply, sewage disposal, 
gas and petrol supply, electrical and telephone communications etc. may suffer 
severe damages due to any disruption caused to the under ground installations 
including buried pipes. 

There are three main causes of failures of buried pipelines due to earthquake; 1) 
excessive axial and bending stresses and deformations in the pipelines created 
mainly due to the phase difference and change of wave shape between different 
stations, 2) fault movement during an earthquake, if the pipeline crosses a major 
fault and 3) inelastic behaviour of soil such as liquefaction. The seismic 
behaviour of buried pipelines is quite different from that of above-ground 
structures (super structures) mainly because a) the inertia forces which are the 
main factors affecting the design of super structures are mainly resisted by the 
surrounding soil in the case of buried pipelines, b) the ground motion is 
considered as coherent for most of super structures, while for buried pipelines 
(being a long continuous system) it is considered as incoherent because of 
phase difference between different stations and change in shape due to the 



variation in soil properties along the pipeline and c) the damage of one super 
structure is restricted to that structure alone in many cases while the damage at a 
certain location within the network of pipes will affect the other portions of the 
system. 


1.2 Literature review 

The available literature on the buried pipes has been reviewed under the 
following categories 

Pipes subjected to static loads 

Pipes subjected to dynamic loads 

1.2.1 Pipes subjected to static loads 

Most of the analytical and numerical methods available so far represent the soil 
system as a continuum. Solution techniques for the problem of soil-conduit 
system can be classified either as analytical methods using classical elasticity 
theory and shell theory to obtain ‘exact’ solutions or numerical methods using 
approximate techniques such as finite difference, finite element and boundary 
element. 

Hoeg (1968) obtained closed form solution of the problem of medium-cylinder 
interaction analytically by developing static solution for an elastic thin cylinder 
embedded in a homogeneous, isotropic and elastic medium with surface loading. 
No-slip condition and full-slip condition at cylinder-medium interface was 
considered. In the solutions for no-slip case, the radial and tangential 



displacements of the cylinder were equated to the corresponding displacements 
of the cylinder at the medium-cylinder interface. In the solutions for the full-slip 
case, the radial displacements of the cylinder and medium were equated, but the 
tangential shear stress was set to zero at the medium-cylinder interface. 

Brown (1967) used the finite element method to analyse the forces on rigid 
culverts under high fills. He used incremental fill sequence technique for the 
analysis, but the soil medium was assumed to be elastic. Moreover, the interface 
boundary condition was treated only approximately by using triangular elements 
which were assigned very high stiffness values compared to fill. Brown et al. 
(1968) improved the representation of the culvert by using three overlapping 
triangular elements. Both rigid and flexible culverts were investigated to study the 
effects of slippage and no slippage. 

Anand (1974) studied the stress distribution around shallow buried rigid pipes 
using finite element method. The study was limited to elastic analysis and depth 
of embedment equal to one diameter of the pipe. The effect of width of imposed 
load, modulus ratio etc. on responses are analysed and compared with closed 
form solutions. Abel et al. (1973) have studied the effect of pipe stiffness, depth 
of cover and effect of slippage at the interface for an elliptical pipe for surcharge 
loading. It was found that the effect of Poisson’s ratio of pipe material and the soil 
is negligible and the slip at the interface is beneficial for establishing arching 
action, particularly for deeply buried pipes. 

Prakash et al. (1976) investigated the buried pipe under embankment using finite 
element method. The parametric studies for height of fill-to-diameter ratio. 



4 


diam 0 ter-to-thickn©ss ratio and backfill propertiss wer© carried out. The effect of 
sequential construction by assuming linear and non-linear behaviour on 
responses was studied. Vailiappan et al. (1977) analysed the buried pipe under 
surface pressure by assuming the medium to be elasto-plastic, limited tension 
material. Conditions of slippage and no-slippage cases are studied by using 
interface elements. Two different yield criteria, namely, Von-Mises and Drucker’s 
modified Von-Mises, have been adopted for the elasto-plastic idealisation of the 
medium. 

Ramakrishnan (1979) analysed buried pipes in linear and non-linear elastic 
media by using finite element method. Pipe-soil system was analysed as a soil- 
structure interaction problem. The pipe-interface-soil was taken as a single unit, 
where, the pipe is simulated by bar elements, interface by interface elements and 
soil by 4-CST quadrilateral elements. Non-linear analysis has been carried out 
simulating stress-strain behaviour by hyperbolic law. The responses like 
springline thrust, crown pressure, crown deflection, springline deflection, 
maximum radial pressure and maximum tangential pressure have been studied 
for various values of parameters, such as radius-to-thickness ratio, modular ratio, 
semi -load width-to-radius ratio for surface pressures. 

1.2.2 Pipes subjected to dynamic loads 

In general depending on the method of analysis and the type of surrounding soil 
the existing literature for dynamic response of buried pipelines, can be classified 
in the following categories 


Damage analysis 



5 


Quasi-static methods 
Analytical methods 
Numerical methods 

1.2.2 .1 Damage analysis 

1.2.2. 1.1 Qualitative damage analysis 

Most of the existing literature concerning buried pipeline damage due to 
earthquake present a qualitative rather than quantitative description of the 
damage. Kachadoorian (1976) reviewed damage data from the 1964 Alaska 
earthquake, the 1968 Mechering earthquake and the 1971 San Fernando 
earthquake and concluded that the geological environment under the buried 
pipeline influenced the intensity and frequency of the pipeline damage. Kube 
(1974), Nasu, Kazama, Morioka and Tamura (1974) have indicated that the 
pipelines moved closely with the ground in both longitudinal and lateral directions 
during the seismic wave propagation. Katayama, Kubo and Sato (1975) 
correlated pipeline damage to pipe size and concluded that smaller pipes are 
liable to break. Luscher, Black and Nair (1975) discussed briefly the design of the 
buried portion of Trans-Alaska pipeline. Duke and Moran (1975) and Whitman, 
Cornell and Tabel-Agha (1975) suggested the use of a reliability/damage level 
approach to pipeline systems design for various intensities of ground motion. 

Ghashi (1977) observed that the damage was maximum in the regions of 
transition from one type of soil to another. Yeh (1977) determined the range of 
conservatism of the assumption that the pipe moves with the surrounding soil. 



The design formulae for pipe stresses due to oblique incident compressional, 
shear and Rayleigh surface waves were derived. 

1. 2.2.1. 2 Quantitative damage analysis: 

In addition to the aforementioned studies, some researchers have analysed the 
quantitative damage of pipelines. Katayma et ai.(1975) (1977) and Katayma 
(1976) analysed the pipeline damage to 1923 Kanto earthquake in terms of 
damage index (D. I). Kubo (1975) studied the pipeline damage from the 19 71 
San Fernando earthquake by dividing the affected regions in to a series of strips. 
He plotted a graph between the number of breaks per unit length versus 
maximum acceleration of ground. Katayma et al. (1977) presented a statistical 
method to analyse seismic damage to 1923 Kanto earthquake employing 
multivariate regression with observed degree of damage. 

1. 2.2.2 Quasi-static method of analysis: 

Sakurai and Takahashi (1969), Newmark (1971) and Shah and Chu (1974) 
analysed the seismic response of underground pipelines to earthquake using the 
quasi-static method without taking into account soil-pipeline interaction, in which 
they assumed that the pipeline was flexible enough to follow the deformation of 
the surrounding soil. Thus by using the maximum values of amplitudes and 
wavelength of the soil seismic deformation, the maximum pipeline strains are 
determined. Nelson and Weidlinger (1977, 1978) suggested the use of an 
‘Interface Response Spectra’ while Wang and Cheng (1978) proposed the use of 
quasi-static approach with soil-pipeline interaction for analysis of buried pipelines 
subjected to earthquake excitation. 
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1. 2.2.3 Analytical methods 

Muleski, Ariman and Aumen (1978) and Novak and Hindy (1978) studied the 
response of buried pipes using shell model. The latter formulated the problem in 
terms of both deterministic and random vibration. Wang and Fung (1979) 
proposed a modified Von-Mises failure criterion to evaluate the safety of buried 
non-homogeneous pipes against earthquakes. Datta et al. (1984) studied the 
seismic response of large cylindrical shells surrounded by an infinite 
homogeneous elastic medium. Two different models 1) a beam model 2) 
Flugge’s shell models are used. In their work, full dynamic interaction between 
the shell and the surrounding soil was taken into account. 

Datta and Mashaly (1986) obtained the response of buried pipelines to random 
earthquake forces using lumped mass model. The earthquake was considered 
as a stationary random process characterized by a power spectral density 
function (PSDF). The cross spectral density function between two random inputs 
along the pipe length was defined in terms of local earthquake PSDF. The soil 
resistance to dynamic excitation along the pipe length has been obtained in 
frequency independent impedance functions which were derived from half-space 
analysis and Mindlin’s static stresses with in the soil due to point loads. In their 
study they considered the cross terms in the soil stiffness and damping matrices. 
Parametric studies have been done for the influence of cross terms in the soil 
stiffness and damping on the response of the pipe. 
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Zerva et. al. (1988) developed a model for near source ground motions, in which 
the excitation was modelled as a random process based on the earthquake 
magnitude and the site characteristics. Power and cross spectral densities of 
acceleration at any station on the ground surface, cross correlation coefficients 
and spectral densities of differential accelerations are the out puts from the 
model. The results of the application of model to the earthquake of January 29, 
1981 in Lotung, Taiwan, has been compared with actual earthquake data records 
at the Strong Motion Array in Taiwan (SMART). The motions obtained by the 
above analytical model are used as input motions for the study of the dynamic 
response of the pipelines. The response of the two different pipeline structural 
models, 1) pipelines with soft joints and 2) continuous pipelines to perfectly and 
partially correlated random input motions in the axial, lateral and vertical 
directions are investigated and the significance of the spatial variation of ground 
motion was also examined. 

Datta and Mashaly (1990) proposed a method based on spatial discretization of 
the pipe with nodal lumped masses for the transverse response of suspended 
spans of submarine pipelines to random excitation. The response was obtained 
in frequency domain using spectral analysis. They have taken into account the 
pressure drag effect and spatial correlation of seismic excitations along the pipe 
length in their parametric studies. 

1. 2.2.4 Numerical methods: 

Numerical methods such as Finite difference method. Finite element method. 
Boundary element method and Hybrid methods were employed by several 



investigators to obtain solutions to dynamic soil-pipe interaction problems. Ang 
and Nev/mark (1963) computed response of the underground structures using 
finite difference method. Farhoomand and Wilson (1970) proposed the non-linear 
finite element method for analysing the blast response of the underground 
structures. Haldor et al. (1980) evaluated the seismic response of submarine 
concrete coated steel pipes buried in a porous bed, using a plane strain finite 
element model. The non linear soil behaviour was treated by an equivalent linear 
method and the effect of energy dissipation in the porous soil medium was taken 
into account by standard viscous boundaries developed by Lysmer and 
Kuhlemeyen (1969). 

Datta et.al.(1984) obtained the displacements and stresses in buried pipeline 
under dynamic loading. A hybrid technique, combining a finite element method 
with the eigen function expressions was used for the dynamic response of pipe. 
The line was modelled as a circular cylindrical shell of small thickness, and the 
incidence disturbances are assumed as plane waves, moving perpendicular to 
the axis of pipeline. Two problems 1) the pipe surrounded by a homogeneous 
soft soil 2) the pipe lying in a cylinder of soft soil which is surrounded by a hard 
material have been studied. 

Yaw-Huei Yeh et.al. (1985) proposed a simplified pipe model with beam on time 
dependent elastic foundation for the dynamic response of buried pipelines in soil 
liquefaction environment. A finite difference method involving both temporal and 
spatial variables was used for the response during earthquakes. The influence of 
a) soil stiffness, b) loading speed of pip6 end, c) initial crookedness and d) 
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supporting length of the pipe segment inside the liquefied region on the pipe 
response are studied. 

Shakib and Abdolsalami (1995) investigated the response of the buried pipeline 
to seismic excitation in lateral and longitudinal directions in time domain. The 
pipeline was idealized with assemblage of 2 D beam element, and a discrete 
lumped mass model was employed to the pipe motions. The soil stiffness and 
damping matrices are derived from static and dynamic continuum solutions. The 
main features of this method are that it can take into account any boundary 
conditions that needs to be imposed at the pipeline ends and also, it can 
consider the cross terms in the soil stiffness and damping matrices. Numerical 
studies are conducted to investigate the effect of pipeline length, element length, 
end conditions of pipe, embedment depth and cross terms of soil stiffness and 
damping matrices in homogeneous medium on the response. 

Stamos and Beskos (1995) have proposed a numerical method (boundary 
element method) for determining the dynamic response of three dimensional 
large underground elastic structures. For transient dynamic disturbance 
numerical inversion of Laplace transform was used for computing response in 
time domain, whereas for harmonic disturbances the response was obtained 
directly in the frequency domain. 

1 .3 Scope of present investigation 

Based on the brief survey carried out, it can be pointed out that an accurate and 
feasible model for seismic response of buried pipeline is not available. The 2 D 
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domain discretization models provide approximate solutions with different levels 
of accuracy. Some of 2 D models with absorbing boundaries provide an 
improvement over the other methods but are not associated with the three 
dimensional earthquake loading. Little research has been done on seismic 
analysis of pipes under general earthquake loading and liquefaction conditions. 

To achieve above objectives, response of buried circular pipes subjected to 
three dimensional seismic loading has been studied using three dimensional 
finite element method in the present investigation. Finite element method based 
package structural analysis programs (SAP-80) has been used in the present 
study. Three dimensional solid elements are used for both discretizing the pipe 
and the surrounding soil. Large three dimensional model, reflect small amount of 
energy through the boundaries has been considered. Modal superposition 
method has been used to evaluate the response of pipe under three dimensional 
earthquake loading. 

1.4 Organization of the Thesis 

The work that has been carried out in this thesis is presented in the following 
manner. 

Chapter 2 describes the general principles and basic concepts of finite element 
method. 3 D descritization and application of finite element method to dynamic 
analysis of pipe-soil system are described. The details of package SAP-80, and 
mode superposition method have been presented briefly. 
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In chapter 3, the results and discussions are presented. Two problems 1) 
response of buried pipe under static loading and 2) response of buried pipe 
under seismic loading have been analysed. The effect of parameters like load 
width-to-radius ratio, radius-to-thickness ratio, depth of embedment-to-diameter 
ratio on responses like horizontal springline displacement, vertical crown 
displacement are studied in case of static analysis. In seismic study, the effect of 
parameters like load width-to-radius ratio, radius-to-thickness ratio, depth of 
embedment-to-diameter ratio, mode shape number, 3 D seismic loading and 
liquefaction zone on responses like horizontal springline displacement, vertical 
crown displacement are studied. 

Conclusions are summarized in Chapter 4. Possible extensions of the present 


work are indicated. 
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Chapter 2 

Formulation and Method of solution 

2.1 General 

The seismic soil-pipe interaction deals with the response of pipe under seismic 
loading. Since soil exists as semi-infinite half space, numerical methods have to 
incorporate the notion of infinity in the formulation. Numerical methods like 
Integral equation method, boundary element method can handle infinite domains 
naturally, where as domain discretization methods like finite difference, finite 
element methods can not be applied to semi-infinite domains directly. A way-out 
of handling such infinite domains is to consider a finite (computational) domain 
for discretization with an approximate boundary conditions (energy absorbing) on 
the finite domain, or a large domain with a boundary which reflects little energy. 

2.2 Finite element method: 

Many engineering problems encountered in the present practice deal with 
systems with distributed parameters. The application of mechanics to such a 
systems often lead to a system of partial differential equations. In case the 
domain and the equations are simple, the solution can be obtained in exact form. 



14 


using Fourier, Laplace transform methods etc. Exact solutions to general partial 
differential equations are difficult to obtain due to irregular and geometrically 
complicated domains. Consequently, various methods of finding suitable 
solutions have been under continuous development. Finite difference method 
and variational methods are more popular ones. The difficulty in applying these 
methods lies in constructing approximating functions of the dependent 
variable, which need to satisfy the geometric boundary conditions on irregular 
domains. The finite element method which is an offshoot of the classical 
methods overcomes these difficulties and gained wide popularity and 
acceptance. 

In finite element method (FEM) the given domain is subdivided in to a finite 
number of small sub domains. The subdomains, called finite elements, are of 
geometrically simple shapes, and permit a systematic construction of the 
approximating functions. Generally these approximating functions are algebraic 
polynomials developed using interpolation theory, and are called interpolation 
functions. Since the approximating functions are defined element-wise, the 
accuracy of the approximation can be improved by increasing the number of 
elements. Generally the interpolation points, called nodes are placed on the 
boundary of the element such that they uniquely define the element geometry. 
Additional nodes which may be required to define the interpolation functions can 
be placed either inside or on the boundary of the element. The boundary nodes 
also facilitate in connecting the adjacent elements together by ensuring the 
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primary degree of freedom to be the same at the common nodes. The nodal 
values of the element displacement uniquely defines the function values within 
the element and they form the basic unknowns of the problem. 

2.3 Formulation 

In this section the basic formulation of finite element method Is discussed. A 
displacement type formulation uses nodal displacements u, v, w as primary 
variables. Variational formulation has been utilized to develop the procedure for 
spatial discretization. Hamilton’s principle has been used to construct the 
functional. The principle is stated as follows (Desai and Abel 1972), 

Among all possible time histories of displacement configurations which satisfy 
compatibility conditions and which also satisfy conditions at time t^ and X 2 , the 
history which is the actual solution makes the Lagrangian functional minimum. 
This can be stated as 

h 

S jldt =0 (2.1) 

'1 

where 6 indicates the variation of the integral and L is the Lagrangian which can 
be expressed as 

L = T-P-W (2.2) 

where T, P and W are kinetic energy, strain energy and potential energy of the 
applied loads. Accordingly 
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L 



p u - e^crj dv + fdv + ^u^rds 
o nr 


(2.3) 


where p denotes the mass density of the material; dots denote the time 
derivative of the variables. 

Lagrangian L of the whole domain can be assembled as sum of the contributions 
over each of the elements, i.e. 




r'' 


(2.4) 


Substituting for the approximate displacement vector w" defined over each of the 
element and stress strain relations in the above equation (Eq. 2.4) results in the 
following expression 




\dv 


(2.5) 


l-{@P 03 ft 0 i- t |-|0 summation symbol is omitted for the sake of conciseness and 
wherever elemental quantities appear, usual summation (assembly) over the 
elements is implied. Substituting the above expressions for Lagrangian in Eq. 2.1 
gives the expression 
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\p{Su^} 


|[A']''{/-}<iv- |[/^]'{r')* 


(2.6) 


di = 0 


Integrating the first term results in the Eq. 2.7, 


I M 

J 


p{Suf 


a 


dt = 


p{Su^] 


t2 


(2.7) 




n 


dt 


Substituting this expression in the Eq. 2.6 results in the following expression 


(2 


n 




dt = 0 


( 2 . 8 ) 


The first term in the right hand side vanishes because the displacement 
configurations should satisfy the boundary conditions at times t^ and t 2 

(Hamilton’s principle) i.e. the variation {Su^} should be equal to zero at times t^ 
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and t 2 rest of the expression inside the parenthesis should vanish because the 
variation is arbitrary. The resulting equation can be written as 

= (2.9) 


where M'l is element mass matrix and is given by 


[m‘]= \p[N^ [N]dv (2.10) 

0 ‘’ 

Assembling the element matrices results in a system of differential equations 
(equations of motion) governing the behaviour the body as follows 

= {/>} (2.11) 

Here the vectors U, and P are time dependent. 


2.4 Finite element model 

An eight noded brick (solid) element is considered as the basic element. Both the 
pipe and the surrounding soil has been discretized with the same eight noded 
elastic elements. The governing equation the motion in which matrices [M] and 
[K] are assembled from the individual element matrices. The element mass and 
stiffness matrices are obtained as described in the following section. 


2.4.1 Element stiffness matrix 

An eight noded isoparametric element is shown in global and local coordinate 
system in Figure 3.1. Since the local coordinates range from -1 to +1, the 
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element is a cube with a side of 2 units long and origin of the system at its 
centroid. The actual element could be any shape with six-faces. The coordinate 
mapping which relates the two coordinate systems can be expressed as 
(Equation 2.13) 




(a) Global 


(b) Local 


Figure 2.1; 8 noded Hexahedron Element 
in Global and Local Coordinate Systems 


x = [7V,]{x,} (2.12) 

where x and Xj are the vectors of independent variables in global and local 
coordinate system respectively, and coordinates for node i respectively; [NJ is a 
diaaonal sub matrix of order equal to the number of space dimensions (3). Each 
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of its entry equals to mapping function of the node. The subscript i varies from 1 
to 8. 

X = [xyzY (2.13) 

' n , 0 o' 

[A^,.]=0iV, 0 (2.14) 

_0 0 N,_ 

=[-T .V, 2,]"' (2.15) 

The mapping function can be expressed as follows (Desai and Abel 1972) 

w, = (1 + , f)(i + n, ’?)(! + i) ( 2 . 1 6) 

where 7, and are coordinates of the node /and take either +1 or -1. 

The displacement vector u can be expressed in terms of its nodal values using 
the shape function expressed in the local coordinate system of the element. 
Since isoparametric formulation is being adopted the shape functions are same 
as the coordinate mapping function of the element. Hence, 

“ = KK“,1 (2-17) 

Using elastic constitutive relations, stresses can be expresses as 

cr =[£■]€ and (T®=[£']e“ (2.18) 

where [E] denotes the elasticity matrix and for three dimensional solid element it 


can be expressed as. 
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Assuming small strain theory, strains can be defined in terms of displacements 
using differential operator [L] as 

& = [L]u (2.20) 

for three dimensional problems the differential operator [L] is given by, 



dx 

0 


0 0 


0 


d_ 
dy 

0 0 — 
dz 

d d 


0 


dy dx 

0 ^ ^ 
dz dy 

^ 0 — 
dz dx 


( 2 . 21 ) 
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Similar relations hold for elemental quantities also, element strains can be 
expressed as 


(2.22) 

where matrix [B] relates element strains to nodal displacement variables, given 
by 

[fi] = [iIW] (2.23) 

The differential operator [L] contains derivatives with respect to physical 
coordinates. The derivatives have to be computed using Jacobian transformation 
from local coordinate space to global coordinate space. The Jacobian can be 
expressed as 













(2.24) 


The ‘, ‘ indicates differentiation with respect to the variable following symbol and 
the symbol Z indicates summation over the number of nodes of the element (8 
here), viz. 


dN, 

dC 


X; + ■ 


dN, 

dC 


•Xj + 


d^ 

dC 


Xj + ■ 


dN, 

dC 


+ 


dN, 


dC 




dN, dN, dN, 


dC 


4" 


dC 


X7 + 


dC 




(2.25) 
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The stiffness matrix of the element now can be expressed as 

[if"] =111 [BWElBl4d(dr,d( (2.26) 

..l-l-l 

where HJH denotes the determinant of Jacobian. 

2.4.2 Element mass matrix 

The element mass matrix can be expressed as (Desai and Abel 1972) 

[m']= j||^,[A'f[W]|./||</fd7d^ (2.27) 

The mass matrix defined by the above Eq. 2.27 is called consistent mass matrix 
because the same shape functions [N] are used to generate the stiffness matrix. 
Often a diagonal form of the mass matrix is adopted assuming the mass of the 
body is concentrated only at the nodes of the body. This form of the matrix is 
diagonal and is called lumped mass matrix. 

2.5 Package used and method of solution 

In this section a brief description of the package used in this study and the 
method of solution has been presented. 

2.5.1 About package 

2.5.1. 1 Structural Analysis Program (SAP-80, 1984) 

It is a general purpose finite element structural analysis program for the linear 
static and dynamic analysis. Its development started as early as 1964 and it was 
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probably or>e of the very first codes based on finite elements. This is the result of 
over ten years of research and development by E. L; Wilson*' at the university of 
California, Berkeley. This program has the capability to analyze very large three 
dimensional systems. With SAP-80 analytical capabilities in terms of elasticity, 
temperature and dynamic analysis are classified into different options. The 
elements available in SAP-80 can be mixed with others within the same option 
analysis Material properties can be (for two dimensional) isotropic, orthotropic, 
temperature dependent, and anisotropic. Boundary conditions can only be 
defined with respect to global Cartesian axes. The input results can be checked 
by plotting. Output results can also be plotted. (SAP-80 manual, 1984). 

2.5.2 Equations of motion 

In earthquake-prone regions, the principal problem of dynamics is the betiaviour 
of structure subjected to earthquake-induced motion of the ground. If the 
displacement of the ground is denoted by u^, the total (or absolute) 

displacement of the mass by u ' , and the relative displacement between mass 
and ground by u (Fig. 3.2), then at each instant of time these displacements are 


related by 

«'(/) = «(/) + », (/) 


(2.28) 
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Figure 2.2 Spring, Damping and Mass System 

Considering the equilibrium of system shown in Fig 3.2 

mil' = -k[ii' (2.29) 

Using the Equations 2.28 and 2.29 and assembling element mass, damping and 
stiffness matrices, the equations of the motion of the system can be obtained as 
M 0 (/) + (.’()(/) + A" U{i) = M ?/^,(/) (2.30) 

where M, C and K are mass, damping and stiffness matrices, respectively. The 
three dimensional relative displacement, velocity and accelerations are indicated 
by U, U and 0 . 

2.5.3 Mode superposition method (Wilson etal. 1981) 

To solve the equation (2.30) mode superposition method has been used in the 
present study. The mode superposition solution involves the introduction of the 
transformation 


t7 = cl>y 


(2.31) 
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where O is the matrix containing three dimensional mode shapes of the system 
and Yis the vector of normal co-ordinates. The introduction of this transform and 


the premultiplication of Equation 2.30 by yields 

+ (2.32) 

For proportional damping the mode shapes have the following properties: 

(f>.^ (2.33) 

K(j)^=cor m.^ ( 2 . 34 ) 

(2.35) 


in which is the /th column of O representing the /th mode shape, m,. is the /th 
modal mass , and 4^, is the damping ratio for mode /. Due to the orthogonality 
properties of the mode shapes, all modal coupling terms of the form A(f)^ are 
zero for / j. Thus, Equation 2.32 reduces to a set of uncoupled equations in 
which the typical ‘modal’ equation is of the following form 

(2.36) 

where 


m. 

is the participation factor for mode /. 

The following equation can be solved for the response y^{t) 


(2.37) 




(2.38) 
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At the point in time where |_v,(/)| is maximum, the response is defined as A 
plot of this maximum displacement versus the frequency for each is by 
definition the displacement response spectrum (wilson et al. 1981) for the 
earthquake Mg(r). 

If the dynamic loading on the structure is specified in terms of displacement 
spectrum, then the maximum response of each mode is given by 

^',ma\ ~ A i.max (2.39) 

Therefore, the maximum contribution of mode / to the total response of the 
structure is 

(2.40) 

For all modes is, by definition, positive. The maximum modal displacement 
L/,- max is proportional to the mode shape , and the sign of the proportionality 

constant is given by the sign of the modal participation factor. Therefore, each 
maximum modal displacement has a unique sign. Also, the maximum internal 
modal forces, which are consistently evaluated from the maximum modal 
displacements, have unique sign. 

In 3D seismic analysis, for combining modal maxima the square-root-of-sum-of- 
squares mehtod causes significant errors (Wilson et al. 1981). A complete 
quadratic combination (CQC) method which reduces errors is used for combining 
the modal responses (Wilson et al. 1981). The CQC method requires that all 
modal response terms be combined by the application of the following equations: 
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For a typical displacement component, u^\ 

(2.41) 

and for a typical force component, : 

A = (2.42) 

where is a typical component of the modal displacement response vector, 
and /i,. is a typical force component which is produced by the modal 
displacement vector, t/,. . 

This combination formula is of complete quadratic form which includes all cross- 
modal terms, and these cross-modal terms may have positive or negative values 
depending on whether the corresponding modal response have the same or 
opposite signs. The signs of modal response are, therefore, an important key to 
the accuracy of the CQC method. 

In general the cross-modal coefficients, , are functions of the duration and 

frequency content of the loading and of the modal frequencies and damping 
ratios of the structure. For a constant modal damping, the expression for cross- 
modal coefficient is 


where 



r'^y + 4(^V(l-l-r)^ 

and ^ is modal damping. 


(2.43) 
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Using the above method of solution and SAP-80 package, response of buried 
pipe under different loads are studied. The results obtained are presented in 
Chapter 3. 
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Chapter 3 


Results and Discussions 


3.1 General 

The formulation described in the previous chapter is applied to solve 1) the 
response of buried pipe under static loading and 2) the response of buried pipe 
under earthquake loading. The basic assumptions are as follows 

1) the soil is an infinite, elastic, homogeneous and isotropic medium. 

2) the pipe is circular in shape , long without joints and elastic. 

3) the pipe is straight and horizontal i.e. there is no level difference between the 
two ends of the pipe. 

4) a perfect bond between pipe and soil i.e. during vibration the displacement of 
pipe must be equal that of adjacent soil. 

3.2 Problem 1: Response of buried pipe under static surface 
loading. 

The response of buried pipe-soil system subjected to static surface loading is 
analysed with the present 3D model. Problem definition and a part of the 
descritized domain are shown in Figures 3.1 and 3.2. L and D are length, and 
diameter of pipe respectively; B and b are half width of the model and half width 
of loading; and H and t are the depth of embedment and thickness of pipe 
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respectively. The geometry and material properties of the problem analyzed are 
as follows 


Size of the model : 7.3 X 10 X 4.8 m 
Number of elements: 256 

Young’s Modulus of the soil medium £^=1X10'^ kPa 
Poisson’s ratio of the soil medium = 0.3 
Young’s Modulus of the pipe £^ = 4.0 X 10® kPa 
Poisson’s ratio of the pipe = 0.2 

Length of the pipe L = 10 m 
Diameter of the pipe D = 0.96 m 
Depth of embedment H = 0.96 or 1 .92 m 

The results obtained are non-dimensionalised as follows, 

the crown (Fig. 3.1) deflections and springline (Fig. 3.1) deflections are non 
dimensionalised by pxR/ where p is the intensity of the surface loading, 

is confined modulus equal to £’^(l-t>j)/(l + L>j)(l-2L>j), and R is the radius of 
the pipe. 

3.2.1 Validation of the model 

The effect of some parameters on the responses like crown displacement and 
springline displacement are studied and compared with the available results 
reported By Ramakrishnan (1979) for validation of the present model. The effect 
of following parameters on responses are studied 

a) variation of responses with load width to radius ratio (b/R) 

b) variation of responses with radius to thickness ratio ( R/t) 
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3.2.1. 1 Effect of load width-to-radius ratio 

The responses of pipe for various b/R ratios have been studied. They are 
studied for a R/t ratio of 20 and Ep/Es of 400. H/D ratios 1 and 2 have been 
considered. 

Crown deflection; The dimensionless crown deflection variation with b/R is 
shown in Fig. 3.3 and Fig. 3.4 . Crown deflection increases with b/R ratio up to a 
value 4, beyond that the deflections tends to be constant. 

Springline deflection; The dimensionless springline deflection variation with b/R 
is shown in Fig. 3.5 and Fig. 3.6 . Springline deflection increases with b/R ratio 
up to a value 4, beyond that the deflections tends to be decreased. 

3.2.1. 2 Effect of radius-to-thickness ratio 

The responses have been studied for various R/t ratios for b/R =1 . The range of 
R/t used are 1 5 to 40. Modular ratio taken for this study is 400. 

Crown deflection; The maximum diametrical deflection occurs at the pipe crown 
relative to invert. The variation of dimensionless crown deflections with R/t ratio 
has been plotted in Fig. 3.7. The crown deflection increases with increasing R/t 
ratio, but the large increase in the thickness changes the crown deflection only 
by a small amount. 

3.2.2 Summary 

Response of buried circular pipe under static loading has been investigated. The 
parameters that effect crown and springline response such as load width-to- 
radius ratio, radius-to-thickness ratio and embedment depth-to-diameter ratio 
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have been studied. Sonne of the important observations made from the study are 
summarized below. 

1. Load width-to-radius ratio has significant effect on both the crown and 
springline response up to a value of 4. Beyond that the crown response 
remains constant, whereas springline response decreases. The maximum 
crown response is about 3.5 times that of the maximum springline response. 

2. The thickness of pipe does not significantly effect the crown response. 

3. Embedment depth has significant effect on both the crown and springline 
responses. With increase of depth of embedment, the response decreases. 
The maximum crown response for H/D = 1 is about 1.3 times that of the 
maximum crown response of H/D = 2. In case of springline response, the 
maximum horizontal springline deflection for H/D = 1 is about 1 .2 times that of 
maximum horizontal springline deflection of H/D = 2. 

The results obtained with the present model are compared with those reported 
by Ramakrishnan (1979) and are shown in Figures 3.3 to 3.7. The results 
obtained compare well with those reported by Ramakrishnan (1979). 




Dimensionless crown deflection, yMs/pR Dimensjonless crown deflection, yMs/pR 
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Fig. 3.3 Dimensionless Crown deflection vs Load width-to-radius ratio 

for H/D = 1 



Fig. 3.4 Dimensionless Crown deflection vs Load width-to-radius ratio 

for H/D = 2 











Dimensionless crown deflection, yMs/pR 
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Fig. 3.7 Dimensionless Crown deflection vs Radius-to-thickness ratio. 
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3.3 Problem 2: Response of buried pipe under seismic loading. 

3.3.1 Validation of the model 

The response of buried pipe-soil system subjected to one dimensional seismic 
loading (in z-direction) is analysed with the present 3D model. Problem definition 
and the descritized domain are shown in Figures 3.8 and 3.9. L and D are 
length, and diameter of pipe respectively; H and t are the depth of embedment 
and thickness of pipe. The geometry and material properties of the problem 
analyzed are as follows- 

Size of the model : 48 X 1 0 X 40 m 
Number of elements; 384 

Young’s Modulus of the soil medium £^ = 8X10'^ kPa 
Poisson’s ratio of the soil medium 0^ = 0.35 
Young’s Modulus of the pipe = 2.0 X 10^ kPa 
Poisson’s ratio of the pipe - 0.2 

Length of the pipe L = 1 0 m 
Diameter of the pipe D = 0.96 m 
Depth of embedment H = 7.5 m 

3.3.1. 1 Input data 


El Centro earthquake (California 1940) data is used as the input and is given at 
the bottom of the model to check the validity of the model. The earthquake data 
and the response spectra of the El Centro earthquake are shown in the Fig. 3.10 
and Fig. 3.1 1 respectively. 



^^0 



Fig. 3.8 Problem definition for 3-D pipe-soil system under seismic loading 



23456789 
Time (sec) 


,11 Response spectrum of 1940 El Centro earthqua 





3. 3. 1.2 Effect of mode shape number 


The effect of mode shape number on the response of crown has been studied, 
and are shown in Fig. 3.12. The vertical crown deflection increases with 
increasing the mode shape number up to a value of 10 and beyond that the 
response remains constant. 

The results obtained are tabulated in Table 3.1 and are compared with the 
results reported by Haidar, A. K. et. al. (1980). 

Table 3.1 Responses under different earthquakes 

Name o f earthquake Intensity of earthquake Vert. Crown response 

Modified 1940 El 

C entro (present model) 0.165g 25.1 mm 

1952 California (Haidar 

model) 0.153g 20.8 mm 





3.3.2 Parametric studies 


The effect of various parameters on the responses like vertical crown 
displacement and horizontal springline displacement are studied with present 3D 
model. The projected descritized model on x-z plane has been shown in Fig. 
3.13. The geometry and material properties used are as follows- 

Size of the model ; 40 X 1 0 X 1 1 m 
Number of elements; 384 

Young’s Modulus of the soil medium £^=1X10"^ kPa or variable 
Poisson’s ratio of the soil medium 0 ^ = 0.3 
Young’s Modulus of the pipe £p = 4.0 X 10® kPa 
Poisson’s ratio of the pipe = 0.2 

Length of the pipe L = 1 0 m 

Diameter of the pipe D = 1 .0 m 

Depth of embedment H = 4.0 m or variable 

The seismic data shown in Fig. 3.10 and Fig. 3.11 are used as input and the 
effect of the following parameters on the responses have been studied. 

a) Effect of three dimensional seismic loading 

b) Effect of mode shape numbers 

c) Effect of radius-to-thickness ratio 

d) Effect of Ep/Es ratio 

e) Effect of embedment depth 

f) Effect of liquefied soil zone 




3.3. 2.1 Effect of three dimensional seismic loading 

The effect of three dimensional loading on the responses are obtained by 
applying the seismic loading (Fig. 3.10, 3.11) in three mutually perpendicular (x, 
y. z) directions to the model shown in Fig. 3.13. The vertical crown displacement 
and the horizontal springline displacement are tabulated in Table 3.2 and are 
compared with the responses of 1-D (applied in z direction only) seismic loading. 
The values of Ep/Eg = 250, R/t = 20, H/D = 4 and V = 1 5 are taken. 


Table 3.2 Effect of 3-D seismic loading on responses 


Crown resp. under 

1-D loading 

mm 

Crown resp. under 

3-D loading 

mm 

Springline resp. 

under 1-D loading 

mm 

Springline resp. 

under 3-D loading 

mm 

19.08 

38.71 

1.82 

24.25 


3. 3. 2. 2 Effect of mode shape number 


The responses of pipe for various mode shape numbers have been studied. 
They are studied for a R/t ratio of 20, Ep/Eg of 250, and H/D ratio of 4. 

Springline displacement: The springline displacement variation with mode 
shape number is shown in Fig. 3.14. The horizontal springline displacement 
remains almost constant with increasing the mode shape number. 

Crown displacement: The crown displacement variation with mode shape 
number is shown in Fig. 3.14. The vertical crown displacement increases with 
mode shape number up to a value 15, beyorid that the displacement tends to be 


constant. 



3.3.2.3 Effect of E^IE^ ratio 

The responses of pipe for various Ep/Eg ratios have been studied. They are 
studied for a R/t ratio of 20, mode shape number V of 1 5, and H/D ratio of 4. 
Springtine displacement: The horizontal springline displacement variation with 
Ep/Eg is shown in Fig. 3.15. With the increase of Ep/Eg i.e. decrease of strength 
in soil, increases the horizontal springline displacement. 

Crown displacement: The vertical crown displacement variation with Ep/Eg is 
shown in Fig. 3.15. similar trend of springline displacement has been observed in 
vertical crown displacement case. 

3.3.2.4 Effect of radius-to-thickness ratio 

The responses of pipe for various R/t ratios have been studied. They are studied 
fora Ep/Eg ratio of 250, mode shape number V of 15, and H/D ratio of 4. 
Springline displacement: The horizontal springline displacement variation with 
R/t is shown in Fig. 3.16. With the increase of R/t i.e. decrease of thickness of 
pipel, increases the horizontal springline displacement, but the large increase in 
the thickness changes the springline displacement only a small amount. 

Crown displacement: The vertical crown displacement variation with R/t is 
shown in Fig. 3.16. similar trend of springline displacement has been observed in 
vertical crown displacement case. 


3.3.2.5 Effect of embedment depth 

The responses of pipe for various H/D ratios have been studied. They are 
studied for a Ep/Eg ratio of 250, V of 1 5, and R/t ratio of 20. 
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Springiine displacement; The springline displacement variation with H/D is 
shown in Fig. 3.17. The horizontal springline displacement remains almost 
constant with increasing the H/D ratio. 

Crown displacement: The crown displacement variation with H/D is shown in 
Fig. 3.17. The vertical crown displacement remains constant with H/D ratio up to 
a value 2, beyond that the crown displacement tends to be increasing. 

3.3.2.6 Effect of liquefied soil zone 

Liquefaction is a state of soil when it becomes unstable. When effective stress in 
soil reduces to zero, the soil loses its shear strength then it behaves like a slurry 
and flows like a liquid. Liquefaction is usually associated with the earthquakes. 
The effect of cyclic process of liquefaction of soil on responses has been studied 
by considering a soil zone of 1.5 times the diameter of the pipe around the 
buried pipe. The Young’s Modulus of this soil Eis, taken as 1/10, and 1/20 of the 
original soil modulus value. The responses under these conditions are tabulate 
in Table 3.3. the values of Ep/E* = 250, R/t = 20, V = 15 and H/D = 4 are taken in 
this study. 


Table 3.3 Effect of liquefied soil zone on responses 



With 

with out 

Response 

liquefaction 

liquefaction 


zone 

E,s=Es/10 Eis=Es/20 

zone 

Springline 

displacement mm 

26.17 

28.19 

24.25 

Crown displacement 

mm 

42.63 

44.68 

38.71 



3.3.3 Summary 


Response of buried pipe under seismic loading has been investigated. The effect 
of 3-D seismic loading, mode shape number, R/t ratio, Ep/Es ratio, embedment 
depth and liquefied soil zone have been studied. A summary of the observations 
made from the parametric study conducted are presented below. 

1 . Three dimensional seismic loading has stronger influence on responses than 
the one dimensional seismic loading. The maximum crown response under 3- 
D loading is almost double that of the maximum crown response under 1-D 
loading. In case of springline response, the maximum horizontal springline 
deflection under 3-D loading is about 14 times the maximum horizontal 
springline deflection under 1-D loading (Ref. Tab. 3.2) 

2. The mode shape number has very little effect on springline response, 
whereas in case of crown response it has significant effect up to a value of 
15 . 

3. With decreasing soil strength i.e., with increasing modular ratio, the pipe 
responses are increasing. The modular ratio has greater effect on crown 
response than that of springline response. 

4. The pipe thickness does not significantly effect the vertical crown and 
horizontal springline response. 

5. The embedment depth has little effect on horizontal springline response, but 
in case of crown response it has significant effect beyond H/D = 2. 

6. Liquefaction zone in soil surrounding the buried pipe has significant effect on 
pipe responses. 












Response in mm. Response in mm 










Chapter 4 


Conclusions and Scope for future study 

4.1 General 

The results obtained in the present study for different class of problems have 
already been discussed in the previous Chapter. The conclusions, possible 
extensions of the study and scope are presented in the following sections. 

4.2 Conclusions 

In this investigation the response of buried pipe under different load conditions 
has been studied using finite element method package SAP-80. Pipe-soil system 
is discretized using 3-D solid elements. Effect of different parameters on pipe 
response are studied using mode super position method. 

From the results obtained for different problems solved in Chapter 3, it can be 
seen that the present 3-D model was more effective in determining the response 
of pipe under general earthquake loading. 

4.3 Scope for future study 

Few possible extensions of the present study could be in the following directions. 

1 . Incorporation of non-linear and non-homogeneous material models. 

2. Incorporation of interface (contact) elements to account for slippage of pipe- 
soil system. 

3. Effect of shape of pipe and elevation difference of pipe ends. 

4. More detailed study of liquefaction effect on pipe response. 
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